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Abstract 

We consider a model Hamiltonian describing the interaction of a single-mode radiation field 
with the atoms of a nonlinear medium, and study the dynamics of entanglement for specific non- 
entangled initial states of interest: namely, those in which the field mode is initially in a Fock 
state, a coherent state, or a photon-added coherent state. The counterparts of near-revivals and 
fractional revivals are shown to be clearly identifiable in the entropy of entanglement. The "overlap 
fidelity" of the system is another such indicator, and its behaviour corroborates that of the entropy 
of entanglement in the vicinity of near-revivals. The expectation values and higher moments of 
suitable quadrature variables are also examined, with reference to possible squeezing and higher- 
order squeezing. 
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I. INTRODUCTION 



A problem of considerable interest in quantum dynamics is that of the identification of 
signatures of non-classical effects in the temporal behaviour of quantum mechanical expec- 
tation values in nonlinear systems. The dynamics of a quantum wave packet governed by a 
nonlinear Hamiltonian provides adequate scope for such an investigation to be carried out, 
as a wide variety of non-classical effects such as revivals and fractional revivals [1], as well 
as squeezing, are displayed by the wave packet as it evolves in time. 

While a generic initial wave packet \ip{0)) governed by a nonlinear Hamiltonian spreads 
rapidly during its evolution, it could return to its original state (apart from an overall phase) 
at multiples of a revival time T^ev, under certain conditions. This is signalled by the return of 
the overlap C{t) = \ {ilj{0)\tlj{t)) p to its initial value of unity ait = nT^-ev Further, at specific 
instants of time in between two successive revivals, fractional revivals of the wave packet 
may occur. This is characterised by the splitting up of the initial wave packet into a number 
of spatially distributed sub-packets, each of which is similar to the original wave packet. 
Both revivals and fractional revivals of a wave packet arise due to very specific quantum 
interference properties between the basis states comprising the original wave packet P]. 

It is evident that, at exact revivals of any initial state, quantum mechanical expecta- 
tion values of observables return to their initial values. Distinctive signatures of different 
fractional revivals show up in the time dependence of the higher moments of appropriate 
operators. By tracking the time evolution of various moments of certain operators, selective 
identification of different fractional revivals can be achieved P]. Since the initial state of 
the system turns out to play a crucial role in determining its subsequent dynamics and the 
non-classical features it exhibits, these signatures also help assess the degree of coherence of 
the initial state '4]. Further, it has been shown that the squeezing and higher-order squeez- 
ing properties of certain quadrature variables in the neighbourhood of a fractional revival 
of the wave packet provide quantifiable measures of the departure from perfect coherence 
of the initial state l5||. The studies in Refs. have been carried out in the context of 

the propagation of a single-mode electromagnetic field in a Kerr-like medium, modelled by 
the Hamiltonian H = hxa^'^O''^, where a , a"^ are photon annihilation and creation operators, 
and X (> 0) represents the susceptibility of the nonlinear medium. (We shall refer to this as 
the case of "single-mode" dynamics in what follows.) The initial state of the field has been 
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taken to be a member of the family of photon-added coherent states as the properties of 
such states include a quantifiable degree of departure from perfect coherence, sub-Poisson 
statistics (a standard deviation that is asymptotically proportional to a power of the mean 
that is less than i), and phase-squeezing. The standard oscillator coherent state (CS) is a 
limiting case of a photon-added coherent state (PACS). 

A similar investigation of the dynamics of two interacting modes is of special interest, 
in view of an additional phenomenon that occurs in this case — namely, entanglement. In- 
teresting aspects of the entanglement that arises when an initial single-mode coherent state 
passes through a nonlinear medium modelled by the above-mentioned Hamiltonian, followed 
by an interaction with a 50% beam splitter, have been discussed by van Enk Sanz et 
al. have examined the non-classical effects that arise in the dynamics of two entangled 
modes governed by a nonlinear Hamiltonian, in the framework of an exactly solvable case: 
two modes of an electromagnetic field interacting in a Kerr-like medium. Taking the ini- 
tial state to be a direct product, either of two Fock states or of two coherent states, the 
periodic exact revival of these states has been established, and the manner in which these 
properties are mirrored in the collapse and revival phenomena displayed by the expectation 
values of appropriate observables has been investigated. The collapses are marked by expec- 
tation values remaining constant over a certain interval, while the revivals are signalled by 
rapid pulsed variations of the expectation values. These features are the close analogues of 
those displayed in the single-mode case mentioned earlier. It must be noted, however, that 
the specific Hamiltonian considered in these studies is symmetric in the two modes, with 
identical nonlinear terms and a symmetric coupling between the modes. As a consequence, 
the Hamiltonian is readily diagonalised, and the resultant dynamics displays a considerable 
degree of regularity, including the occurrence of revivals. 

However, this symmetry is not expected to be present in a generic Hamiltonian governing 
the interaction of a single-mode field with a nonlinear medium. The best one can hope for 
is the occurrence of approximate or near-revivals for most values of the parameters in the 
Hamiltonian. In addition to near-revivals and related phenomena, several other features of 
interest are exhibited in general by the dynamics of two interacting modes. Our objective in 
this paper is to demonstrate and investigate these. We will also examine the link between the 
extent to which revivals occur and the nature of the initial state. A suitable Hamiltonian 
for our purposes is the one -9] that describes the interaction of a single-mode field with 
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the atoms of the nonhnear medium through which it propagates. The latter is modelled 
by an anharmonic oscillator. The important point to note is that this Hamiltonian is not 
symmetric in the two interacting degrees of freedom, and, moreover, is not diagonalisable in 
general. 

In order to bring out the salient features of entanglement dynamics, we consider initial 
states that are direct products of the field and atom modes. Entanglement of these two modes 
sets in during temporal evolution. However, for certain ranges of values of the parameters 
and coupling constants in the governing Hamiltonian, we show that the modes disentangle at 
specific instants during the evolution, and the state of the system returns close to its initial 
form (apart from a phase). This feature is clearly the analogue of wave packet revivals 
in the dynamics of a single mode. Appropriate indicators to determine quantitatively the 
extent and nature of near-revivals and their fractional counterparts are the entanglement 
entropies, as measured by the sub-system von Neumann entropy (SVNE) and the sub- 
system linear entropy (SLE), the sub-system considered here being the field mode. These 
entropies show marked dips at near-revivals, and are also seen to reduce significantly in 
certain cases at fractions |, |, | and ^ of the near-revival time Trev, signalling the appearance 
of counterparts of fractional revivals. We have also examined the behaviour of another 
indicator of near- revivals, the "overlap fidelity" C(Trev), as a function of a relevant parameter 
in the Hamiltonian, namely, the coupling between the sub-systems represented by the two 
modes. This fidelity is defined as the maximum value that the overlap | {■ip{0)\tp(t)) p attains 
in the vicinity of the first near-revival, for a given initial state. We study the role of the 
specific (non-entangled) initial state considered, and the extent of departure of the initial field 
mode from coherence, on the subsequent revival properties of the state. In particular, the 
link between the entanglement of states and their squeezing properties has been investigated. 
The initial state of the field mode is taken to be a member of the family of photon-added 
coherent states. The advantages of a PACS have already been stated. We add that this 
family of states is now within the realm of experimental realisation in the foreseeable future, 
a single-photon added coherent state having recently been generated experimentally and 
characterised using quantum tomography 10]. With increasing departure from coherence 
of the initial field mode, the entropy of entanglement at any instant during the temporal 
evolution of the quantum state of the system also increases, and even near-revivals do not 
occur. 
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The plan of the rest of this paper is as follows: In Section m we discuss the relevant fea- 
tures of the model Hamiltonian (ol we use to study the dynamics of two- mode entanglement. 
In Section IIIIl we examine three different indicators or measures of the extent to which a 
state revives at any instant: the first of these comprises the instantaneous SVNE and SLE. 
The trends observed here are corroborated by the behaviour of the overlap fidelity, which 
we examine as a function of the strength of the coupling between the field and the medium. 
Finally, certain operators whose expectation values carry signatures of the different frac- 
tional revivals are also identified. We have examined the dynamics of several initial states, 
taking the atom to be in the ground state while the field is, respectively, in a Fock state, 
a CS, and a PACS. This enables us to analyse systematically the effects of different initial 
field modes on the dynamics. The relationship between the squeezing property of the state 
of the system and the initial field mode is also brought out. 



II. SINGLE-MODE FIELD IN A NONLINEAR MEDIUM 

The interaction of a single-mode field of frequency uo with the atoms of the nonlinear 
medium through which it propagates is modelled by the HamiltonianQ 

H = uj a^a + u^b'^b + -fb^%'^ + g {a% + b^a) . (1) 

(We have set ^ = 1.) a and are the annihilation and creation operators pertaining to 
the field, while b and 6^ are the corresponding atom operators. The medium is modelled by 
an anharmonic oscillator with frequency uoq and anharmonicity parameter 7. The coupling 
constant (7 is a measure of the strength of the coupling between the field mode and the atom 
mode. It is easily verified that the total number operator Ntot = + b'^b commutes with 
H . We reiterate that the foregoing Hamiltonian is not symmetric in the two modes, even 
when oj = ojq. 

The Fock basis is given by {\n')^ ® I^t-);,}, where n' and n are the eigenvalues of a^a and 
6^6, respectively. In the absence of the coupling constant g, H is trivially a direct sum of 
Hamiltonians that are functions of the number operators for the two modes. In the absence 
of the anharmonicity parameter 7, the coupled Hamiltonian is essentially linear in each of the 
sub-system variables, and can be diagonalised in terms of linear combinations of the original 
ladder operators. In physical terms, this leads to fairly simple dynamics, essentially entailing 
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a simple periodic exchange of energy between the two sub-systems or modes. When both g 
and 7 are non-zero, the system displays a wide variety of dynamical behaviour, depending 
on the value of the ratio 7/(7. Additional insight into the nature of the model Hamiltonian 
is gained by re-expressing it in terms of angular momentum operators defined in the usual 
manner, according to J+ = a^h , J_ = , = |(a^a — We then have 

H = h [liuj + UJO- 7)Ntot + {uj-uJo + 7) J, + i7(Ntot - 2J,f + 2gJ^] , (2) 

where Ntot (Ntot + 2) = 4J^. As [Ntot , H] = 0, we may write the basis states as — n)^ ® 
\n)^, using N to label the eigenvalues of N^t • For notational simplicity, let us write 

\N -n)^(g)\n)f^= \N -n; n) . (3) 

It is evident that (A^ — n; n\H \N' — n' ; n') = 0, if A^ 7^ A^'. Hence, for each given value 
of A^, the Hamiltonian H can be diagonalised in the space of the states \N — n \ n), where 
n = 0, 1, . . . , A^. Let the eigenvalues and eigenstates of H be X^s and \ipNs), respectively, 
where s = 0, 1, . . . , A^ for a given A^, and A^ = 0, 1, ... ad inf. It is convenient to expand 
lipNs) in the basis {\N — n; n)} as 

N 

\M = Y.'^n' \N-n;n), d^' = {N - n; n\^^s). (4) 

n=0 

In this basis, (each block of) H can be represented as a real, symmetric, tridiagonal matrix. 
Its eigenvalues X^pand eigenstates \ipNs) can be found numerically using appropriate matrix 
algebra routines Q^Q]- An initial state {ipiO)) of the system evolves in time to the state 

00 N 

m)) = m \m) = E E i-^^Nst) {^Ns\m) \^ns) (5) 

N=0 s=0 

at time t. For our purposes, it is necessary to compute the time-dependent density operator 
pit) for different choices of the initial state \ip{0)), as well as the reduced density matrices 
corresponding to the field and atom sub-systems. The main steps in the procedure are 
outlined in the Appendix. 

As already mentioned, if either (/ or 7 is equal to zero, the Hamiltonian in Eq. is exactly 
solvable, and there is periodic exchange of energy between the field and atom oscillators. For 
non-zero values of the ratio 7/(7 of the respective strengths of the nonlinearity and the field- 
atom interaction, collapses and near-revivals could occur over certain intervals of time, in 

6 



between these periodic exchanges of energy. This phenomenon translates into the behaviour 

of expectation values of certain observables as well. For instance, during a collapse of the 

energy exchange over an interval of time, the mean photon number {|a^a) remains essentially 

constant. A revival of the energy exchange is signalled by rapid oscillations of the mean 

photon number about this value, over the relevant time interval. 

Further, when the atomic oscillator is initially in its ground state, while the field starts 

either in a Fock state or in a coherent state, one finds the following resultsj^: (a) For 

weak nonlinearity (7/(7 <^ 1), collapses and revivals of the mean photon number occur 

almost periodically in time, for both kinds of initial field states. The near-revival time is 

I I 

approximately given by 27r/7 in the former case, and 47r/7 in the latter [13]. (When 7 
is exactly equal to zero, there is no nonlinearity in H, and the system is merely periodic. 
There is no question of revivals in this case.) (b) For 7/(7 ~ 1, such collapses and revivals 
occur more irregularly if the field is initially in a coherent state, compared to the case when 
it is initially in a Fock state. As the nonlinearity is increased (7/(7 ^ 1), collapses and 
revivals gradually become less discernible. Bearing these results in mind, in the next section 
we examine the manner in which the above-mentioned collapse and revival phenomena are 
mirrored in the entropy of entanglement of the system. We identify suitable observables 
which carry signatures of collapses and revivals, and discuss the influence of the departure 
from coherence of the initial state of the field on the extent to which it revives subsequently. 

III. ENTANGLEMENT PROPERTIES 

We now examine the detailed dynamics of three different initial states which are direct 
products of the field and atom states, evolving under the Hamiltonian in Eq. (^J). As stated 
earlier, the initial state of the atom is taken to be the oscillator ground state \0)^, while that 
of the field is, respectively, (a) a Fock state \n')^ , n' = 0, 1, . . . ; (b) a CS \a)^ ; and (c) an 
m-photon-added CS \a, m)^ , m = 1, 2, . . . . (Recall that the suffixes a and h correspond to 
the electromagnetic field and the atoms of the medium, respectively.) The CS and PACS 
referred to have the standard expansions in the Fock basis, namely. 




(6) 
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and 



a, m) = 



(at)™ 


\a) 


(at)™ 


\a) 




a™ at™ 




(-^) 



(7) 



where a G C , = lap, and is the Laguerre polynomial of order m. 

As we are dealing with a pure bipartite system, it is natural to consider the time- 
dependences of Sk , the sub-system von Neumann entropy (SVNE), and 6k , the sub-system 
linear entropy (SLE), where the suffix k stands for either a or b, depending on the sub-system 
considered. These quantities are defined as 

Sk{t) = -TTk{pk{t)lnpk{t)} (8) 

and 

5,(t) = l-Tr,{p2(t)}, (9) 

where pk{t) is the time-dependent reduced density operator for the sub-system concerned. 
In terms of the set of eigenvalues of Pk(t), we have 

Skit) = A«(t) In A«(t), 6k{t) = 1-Y, {^t\t)}\ (10) 

i i 

where the summation runs over all the eigenvalues X^l\t). 

n 

Our objective is to investigate the detailed dynamics of entangled states jl^ exhibiting 
revival phenomena. We therefore restrict ourselves here to the case of weak nonlinearity, i. e., 
7/(7 << 1, as this is the situation in which these phenomena occur most unambiguously. For 
illustrative purposes, we set the values of the parameters aX u = ujq = 1, and 7 = 1, g = 100, 
so that 7/(7 = 10~^. As stated in the preceding section, some of the steps involved in the 
calculation of the foregoing sub-system entropies are given in the Appendix. For pure states 
in a bipartite system, of course, Sa = Sb and 6a = 6b at any instant of time. Figures [HJa) 
and (b) depict plots of SVNE and SLE versus gt for respective initial states |10)^ <S) |0)^ (or 
|10; 0), in our notation) and \a)^ ® |0){, = |a; 0) with the parameter value i/ = 1. The 
band-like appearance of the plots in Fig. ^a) arises from the extremely rapid oscillations 
of the ordinates. The corresponding plots for an initial state \a,5)^ ® |0)^ = |(a,5) ; 0) in 
the cases u = 1 and u = 5 are shown in Figs. Efa) and (b) respectively. In all the cases 
above, the SVNE (the upper plot in each figure) is larger than the SLE (the lower plot in 
each figure) at any instant of time. It is evident that both SVNE and SLE display roughly 
similar oscillatory behaviour in time. However, certain striking differences arise in the time 
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evolution of the SVNE and SLE, depending on the actual initial state considered. If the 
field is initially in a Fock state or a CS, the entropies return to values close to zero at regular 
intervals of time (see Fig. P), signalling a near- revival of the initial state. (Note that the 
revival times are indeed approximately equal to 2tt and An, respectively, recalling that we 
have set 7 = 1 and g = 100.) In contrast, if the initial state of the field is a PACS, the 
extent of revival is considerably reduced (see Figs. Efa) and (b)). Further, with an increase 
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in the value of u, the oscillations in the SVNE and SLE die down. This effect is enhanced for 
larger values of m, as seen in the rapid increase and saturation of both the SVNE and SLE 
for an initial state \{a, 5) ; 0) for z/ = 5, in contrast to the corresponding plots for u = 1. 

It is also clear that the SVNE and SLE display marked oscillatory behaviour near 
|Trev, l^rev and ^Trev This bchaviour may be regarded as the counterpart, in our coupled 
system, of the fractional revivals seen in the case of a single-mode nonlinear Hamiltonian 
3, Again, these oscillations die down in amplitude with increasing m when the initial 
state of the field is a PACS, and are most pronounced if the field is initially in a Fock state. 

As mentioned in Section 1, another indicator that characterises the degree of revival of 
an initial state is the overlap fidelity, defined as 



C(T,ev)= max | (^/'(O)lV'(t)) 



(11) 



which is the maximum value attained by the overlap function in the vicinity of the first 
near-revival time corresponding to a given initial state \iIj{0)). In Figs. Efa) and (b), we 
have plotted this quantity as a function of g, the strength of the coupling between the field 
and atom modes, again for non-entangled initial states in which the field is in a CS or a 
PACS. We see quantitatively how, with an increase in the departure from coherence of the 
initial field state, near-revivals occur only for ever increasing values of the coupling strength 
g relative to the coefficient 7 of the nonlinearity in the Hamiltonian. 
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FIG. 3: Overlap fidelity vs. g for an initial state |(a,m) ; 0), with 7 = 1 and (a) = 1 and (b) 
1^ = 5. 
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Next, we consider whether signatures of the features seen above appear in the time 
evolution of expectation values of observables. For this purpose, we define the quadratures 



(12) 



where 



and 



Xn = [a 



(13) 



= (a - at)/(zV2), pb = {b-b^)/{tV2). (14) 
If the field is initially in a Fock state, both (^) and (77) vanish identically at all times. Setting 
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FIG. 4: (^) vs. gt for an initial state |a;0) with and u = 1 and u = 5, respectively (7/5 = 10 ^). 



uj = ujq = 1, '-f = 1 and g = 100 as before (for ready comparison with the time evolution 
of the SVNE and SLE discussed above), we have plotted {^) versus gt for an initial state 
\a; 0) with z/ = 1 and 5, respectively, in Figs. IH^a) and (b). Figure IH^a) shows that (^) 
displays rapid pulsed oscillations near t = Air, similar to its behaviour near t = 0. This 
manifestation of revivals is consistent with the behaviour of the SVNE and SLE in this case, 
cf. Fig. ^b). The collapses are not sharp, in the sense that (^) is not constant over the 
time interval between successive revivals — oscillatory bursts occur in between, with a slight 
enhancement of these oscillations around the fractional revival at ^Trev In contrast, the 
collapses in between revivals are much more complete for larger values of u, as seen in Fig. 
Efb), consistent with the corresponding behaviour of the SVNE and SLE in this case. An 
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FIG. 5: (0 vs. gt for an initial state |(q, 1) ; 0) with u = 1 {-j/g = 10 ^). 



interesting feature is that these collapses become much sharper for even a marginal departure 
from coherence of the initial state of the field (i. e., even for as low a value as m = 1). (^) 
remains virtually constant over the duration of the collapse, and then bursts into rapid 
oscillations close to revivals, as seen in Fig. El which corresponds to m = 1 and u = 1. As 
in the case of single-mode dynamics, the amplitude of the oscillations in the neighbourhood 
of Trev decreases significantly with an increase in m. Thus, for small values of z/, it is easy 
to distinguish between an initial CS and an initial PACS. The expectation value of 1] also 
displays these signatures. We may add that, while the sub-system variables Xa, Xb, Pa and pb 
do exhibit near-revivals in their expectation values, their higher moments do not capture the 
occurrence of fractional revivals. However, the higher moments of the combinations ^ and 
rj do carry distinguishing signatures to selectively pin-point the analogues of the different 
fractional revivals that occur in the single-mode case. Hence ^ and t] are the appropriate 
dynamical variables in the interacting system under consideration. 

The standard deviation A,^ of ^ reflects the occurrence of the dips in the plots of the 
SVNE and SLE at ^Tj-cv The plot of versus gt for an initial state |10 ; 0) shows a burst 
of rapid oscillations at t ~ tt (Fig. This feature holds for an initial CS or PACS as well, 
as is evident (Fig. [7j) from the sudden burst of oscillations in A^ around t ~ 27r for initial 
states \a; 0), |(a, 1) ; 0) and |(a, 5) ; 0) (recall that T^^v — 4vr in this case). 

We note that squeezing occurs in the neighbourhood of |Ti.ov when the initial state of the 
field is a coherent state: A^ drops below the value | (the horizontal dotted line in Fig. Ej) 
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in the case when the initial state is |a ; 0), in contrast to what happens for an initial PACS 
\{a,m) ; 0). While this is similar to the squeezing property seen in the case of single-mode 
dynamics, such parallels do not hold in the case of higher-order squeezing. The relevant 
quadrature variables in this case are obvious generalisations of those considered in the case 
of single-mode dynamics and are given by Q| 



ii + at? + W + b^i) 



i« - a^" + 6' - b^'^) 



(15) 



2V2 ' 2iV2 

In contrast to the single-mode example, even for weak nonlinearity (7/(7 = 10^^) and v = 1, 
amplitude- squared squeezing (g = 2) is absent at t = |Trev whether the field is in a Fock 
state, or a CS, or a PACS. 
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FIG. 6: vs. gt for an initial state |10; 0) (7/5 = 10" 



Finally, turning to the higher moments of ^ and 77, we note that all odd moments of ^ 
vanish identically for all t if the initial state is a direct product of Fock states. For small 
values of z/, the higher moments of ^ show distinct signatures at fractional revivals only if m 
is sufficiently large. However, for larger values of z/, such signatures appear even in the case 
of an initial CS (m = 0). In contrast to this, the temporal evolution of the variances and 
higher moments of the subsystem variables x^, Xf,, Pa and pi, do not display these signatures. 
In this sense, the expectation values of sub-system quadrature variables would seem to be 
inappropriate choices for investigating collapse and revival phenomena in the presence of 
entanglement, although, as we have shown, the sub-system entropies are eminently suitable 
indicators in this regard. 
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FIG. 7: vs. gt for an initial state \{a,m) ; 0) with z/ = 1 and m = 0,1 and 5, respectively 



The model Hamiltonian we have used to study two-mode dynamics leads to several other 
interesting features that are manifested in the expectation values of sub-system variables, 
as the parameters in H are varied. In particular, the ergodicity properties of the system, 
in a "phase space" spanned by such expectation values, exhibit a range of behaviour from 
quasi-periodicity to exponential instability — the latter, notwithstanding the fact that the 
classical counterpart of H is an integrable two-freedom Hamiltonian. These results will be 
reported elsewhere. 
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APPENDIX: CALCULATION OF THE DENSITY MATRIX 

We outline here the procedure used for calculating the density matrix elements required 
for the determination of the entropies and related quantities in the main text. 

From Eq. ^ for the state vector of the system at time t, it follows that the time- 



(7/5 = 10-2). 
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dependent density matrix p(t) is 

oo N oo N' 

pW = X15Z5Z5Z exp [-i(Ajvs - Atv's')^] 

7V=0 s=0 N'=0 s'=0 

X {^nMO)) (^(O)I^^v) I^jvs) {i^N's'l . (A16) 

For instance, if the atomic oscillator is initially in the ground state |0)^ and the field is in 
the Fock state \N)^, the density matrix reduces to 

N N 

P{t) = 5^ 5^ exp [-z{Xms - XNs')t] C C li^Ns) {^Ns'\, (A17) 
s=0 s'=0 

in terms of the expansion coefficients c?^** defined in Eq. (0}. (This expression is explicitly 
A^-dependent, as expected.) Hence 

(tpMilpit) liJM'i') = exp [-i{\Mi - \M'i')t] 4" 4^''' ^inSnm'. (A18) 
It is evident that p{t) is effectively an (A^ + l)-dimensional diagonal matrix in this case. 

For an initial state \{a,m); 0) we find, using Eq. ()A16|) . 

_,, oo Af oo N' 



m\Lm(-iy) ^ ^ ^ ^ iN -m)\ (N' - m)\ 

™^ ' N=m s=0 N'=ni s'=0 ^ ' ^ ' 

X exp[-z(A7Vs - Ajv'sO^lc^o'' c?^''' I^ws) (^w's'h (A19) 

where we have used the expansion of the PACS |a;, m) in the Fock basis. The corresponding 
matrix elements of the density matrix are given by 



\JM\ M'\ (--'^^-'^ l^,*\M'-m 



m\Lm{-Jy) {M - m)\ {M' - m)\ 
X expH(AM/-AM';')^]CC- (A20) 



Here, and in the rest of this Appendix, it is understood that contributions from terms of 
the form l/(— n)!, where n is a positive integer, vanish. 

The expectation values and higher moments of the quadrature variables ^{t) and i]{t), 
defined in Eq. ()12|) . can now be obtained numerically, using the above expressions for the 
density matrix and the matrix elements of the operators a and b in the basis {ipMs)- The 
latter are given by 

N' 

{ijNsl a liJN's') = J2^N' - rff rff ^' s^^r,'-i (A21) 

n=0 
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and 

N' 

ii^Nsl b lijN's') = Yl ^-1 dn''' Sn,N'-1 , (A22) 

n=l 

respectively. As these are purely off-diagonal, and p{t) is diagonal for an initial state that 
is a direct product of Fock states, it follows that all the odd moments of ^ and 77 vanish 
identically for all t, as asserted in the text. This is no longer true for the other classes of 
initial states considered. 

The time-dependent reduced density matrices Pk{t) {k = a,b) are given by 

00 

Pa{t) = TT,[p{t)] = J2b{n\p{t)\n),, 

n=0 

00 

p,it) = IV„ [p{t)] = Yl Pi^) l^)a ■ (A23) 

n=0 

Corresponding to an initial state ; 0), these reduced density matrices Pk{t) take the form 

N N N 
n=0 s=0 s'=0 

X C C C C \{N - n))„ „((iV - n)\ (A24) 

and 

N N N 

Pb{^) = ^ exp [-i(AAr5 - Aiv^O^] 

n=0 s=0 s'=0 

X C C <in 4-n I (A^ - rz)), ,((7V - n) I . (A25) 
Hence we have, in the Fock basis, 

N N 

a{n\ Pa{t) = ^ ^ exp [-i{\Ns - >^Ns')t] C C"' ^N-n ^N-n' ^nn' 

s=0 s'=0 

(A26) 

and 

AT N 

M p,{t) |n\ = 5^ 5^ exp [-i(A^, - Ajv.O * ] C C C C <^nn' ■ 

s=0 s'=0 

(A27) 

As before, these are explicitly A^-dependent finite- dimensional matrices. 
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For an initial state \{a,m) ; 0) the expressions for pk{t) are given by 



oo oo oo N 



^ - V V V V V v^^^^ («)^~" (^l^'"'" 

^"^^ m\Lm{-y)^ ^ ^ ^ ^ iN -m)\(N' -m)\ 

X expH(A^,-A^v)t]C<'^'rff \{N -n))^,{iN' -n)\ 



(A28) 



and 



oo oo Af oo N 



^ VWYNn. (ar-"^ {a*r- 

P'^^> - ^\T.(-.\2^ 1^ 1^ 1^ 1^ 



(A29) 

where A^min = max(n, m). These are infinite-dimensional matrices. The corresponding 
matrix elements of pA;(t) in the Fock basis are given by 

e- y !^ v/(n + /)! (n + /Q! 

a^MPaWIWa 2^ 2^ + / _ ^) 1 + // _ ^) I 

^ ' "min S = S' = ^ I \ ' I 

(A30) 

and 



//I lA ^ \- \- \- ^(n + /)!(n + /')! («)"-^'-™ 



m!L^(-z/) ^^^^ (n + Z-m)!(n + Z'-m)! 

X exp [-^(A^. - \^,,)t\ dt"-'^' d^t'^'' dt"^'^' 4"+'')^ 

(A31) 

where rimin = max (m — l,m — I'). The corresponding expressions in the case of an initial 
state |a ; 0) are obtained from the above by simply setting m = 0. 

All the reduced density matrices (with elements given by Eqs. (fX2fi|l . (jMU, |A8n|l and 
()A31|) ) are hermitian. They are diagonalised numerically, and their eigenvalues are used to 
compute the entropies Sk{t) and Sk{t) given by Eqs. (fTUI) . In the case of infinite-dimensional 
matrices, rapid convergence in numerical computation is provided by the factorials in the 
denominators of the summands in the expressions derived above for the matrix elements. 
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We use double precision arithmetic with an accuracy of 1 part in 10^. As mentioned in the 
text, we use the equahty of Sait) and Sb(t), and the condition p^(t) = p(t) for the density 
matrix of the total system, as some of the checks on the numerical computations, as we are 
only dealing with pure states. 
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